Draft version May 12, 2009- VERSION 23 
Preprint typeset using 1^1^^ style cmulatcapj v. 05/04/06 



INDEPENDENT CONFIRMATION AND REFINED PARAMETERS OF THE HOT JUPITER XO-Sb^ 



o 
o 

>>■ 



6 



A. Pal^'-'-'', G. a. Bakos^-"', J. Fernandez-, B. Sipocz^-^ G. Torres^, D. W. Latham^, Geza Kovacs"*, R. W. NoYES^ 
G. W. Marcy", D. A. Fischer', R. P. Butler**, D. D. Sasselov^, G. A. EsQUERDO^ A. Shporer", T. Mazeh'^, 

R. P. Stefanik^, H. Isaacson" 

Draft version May 12, 2009- VERSION 23 

ABSTRACT 

We present HATNet observations of X0-5b, confirming its planetary nature based on evidence 
beyond that described in the announcement of Burke et al. (2008), namely, the lack of significant 
correlation between spectral bisector variations and orbital phase. In addition, using extensive spec- 
troscopic measurements spanning multiple seasons, we investigate the relatively large scatter in the 
spectral line bisectors. We also examine possible blended stellar configurations (hierarchical triples, 
chance alignments) that can mimic the planet signals, and we are able to show that none are consistent 
with the sum of all the data. The analysis of the S activity index shows no significant stellar activity. 
Our results for the planet parameters are consistent with values in Burke et al. (2008), and we refine 
both the stellar and planetary parameters using our data. X0-5b orbits a slightly evolved, late G 
type star with mass A/* ~ 0.88±0.03Mo, radius i?^, = 1.08 ± 0.04 i?©, and metallicity close to solar. 
The planetary mass and radius are 1.059 ± 0.028 Mj and 1.109 ± 0.050 respectively, corresponding 
to a mean density of 0.96^°;" g cm'^. The ephemeris for the orbit is P = 4.187757 ± 0.000011 d, 
E = 2454552.67168 ± 0.00029 (BJD) with transit duration of 0.1307 ± 0.0013 d. By measuring four 
individual transit centers, we found no signs for transit timing variations. The planet X0-5b is notable 
for its anomalously high Safronov number, and has a high surface gravity when compared to other 
transiting exoplanets with similar period. 

Subject headings: planetary systems — stars: individual (XO-5, CSC 02959-00729) techniques: spec- 
troscopic 
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1. INTRODUCTION 

There are numerous dedicated transit searches survey- 
ing the sky for extrasolar planets that periodically transit 
across the face of their host star. Among the wide an- 
gle searches, those presenting discoveries have been TrES 
(Brown & Charbonneau 2000; Dunham et al. 2004), XO 
(McCuUough et al. 2005; Burke et al. 2007), HATNet 
(Bakos et al. 2002, 2004), and Super WASP (PoUacco et 
al. 2006; Cameron et al. 2007). The initial high hope 
of finding hundreds of such planets (Horne 2001) was 
followed by 5 years of poor harvest, and a steep learn- 
ing curve for these, and many other projects. In retro- 
spect we now understand that several important factors 
had initially been underestimated, such as the need for 
dedicated telescope time, optimal precision, stable in- 
strumentation, low systematic noise, the number of false 
positives (Brown 2003), optimal follow-up strategy, and 
access to high precision spectroscopic instruments. The 
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last year showed an exponential rise in announcements^", 
indicating that these dedicated efforts have started to 
bear fruit. In fact, they have reached a success rate such 
that the same object is occasionally independently found 
and announced by different groups (WASP-llb: West et 
al. (2008) = HAT-P-lOb: Bakos et al. (2008)). Such sce- 
narios are not necessarily duplication of effort. It is re- 
assuring that completely independent discoveries, follow- 
up observations and analyses lead to similar parameters. 
They also provide an opportunity for joint analysis of all 
datasets. Here we report on a similar case, the confirma- 
tion of the planetary nature of the transiting object XO- 
5b, announced by Burke et al. (2008). The present pa- 
per provides not only strong new evidence supporting the 
planetary nature of the object, but also improved phys- 
ical properties that aid in the comparison with theories 
of planet structure and formation. In § 2 we describe the 
details of the photometric detection. The follow-up ob- 
servations, including the discussion of the bisector span 
measurements are presented in § 3. The subsequent steps 
of the analysis in order to characterize the star, orbit and 
the planet are discussed in § 5. 

2. PHOTOMETRIC DETECTION 

Two telescopes of the HATNet project, namely 
HAT-6, stationed at Fred Lawrence Whipple Observa- 
tory (FLWO, A = 111°W), and HAT-9, located on the 
rooftop of the Submillimeter Array control building at 
Mauna Kea, Hawaii (A = 155°W), were used to observe 
HATNet field "G176" (a = 07*^28'", 6 = -l-37°30') on a 
nightly basis between 2004 November 26 and 2005 May 

-'^^ http://www.oklo.org, http://www.exoplanet.eu 
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9. Altogether we acquired with these telescopes 2640 and 
4280 frames, respectively, with exposures of 5 minutes. 

A number of candidates have emerged from this field, 
and have been subjected to intense follow-up by larger in- 
struments (§3). One candidate has become the transit- 
ing planet we call HAT-P-9b (Shporer et al. 2008). An- 
other candidate internally labeled HTR176-002 has re- 
ceived extensive follow-up over the past two years. How- 
ever, the large scatter in the spectral line bisectors, and 
their tentative correlation with orbital phase discouraged 
us from early announcement, and motivated us to pursue 
it further. Subsequently. HTR176-002 was annoimced as 
X0-5b by the XO group in 2007 May (Burke et al. 2008, 
hereafter BOS). Nevertheless, we present here our results 
since they provide independent confirmation and also re- 
fine most of the parameters. 

By chance, XO-5 happens to fall at the edge of field 
"G176" which overlaps with field "G177" (a = OS'^OO™, 

5 = -F37°30'). This field has been observed by the HAT- 
Net telescope HAT-7 and by the WHAT telescope at 
Wise Observatory, Israel (Shporer et al. 2006). Using 
these telescopes we collected 5440 and 1930 frames, re- 
spectively. Altogether we obtained ^ 14290 frames with 
photometric information on XO-5 — an unusually rich 
dataset compared to data available for a typical HATNet 
transit candidate. 

The frames from field "G176" were processed and an- 
alyzed as described e.g. in Bakos et al. (2007). The 
light curves from this field were corrected for trends 
using the method of External Parameter Decorrelation 
(EPD, see Bakos et al. 2009), and the Trend Filter- 
ing Algorithm (TFA; Kovacs et al. 2005). The fight 
curves were then searched for periodic box-like signals 
using the Box Least Squares algorithm of Kovacs et al. 
(2002). We detected a significant dip in the light curve 
of the I « 12.17 magnitude star GSC 02959-00729 (also 
known as 2MASS 07465196+3905404: a = 7'^46'"51^97, 

6 = -F39°05'40".5; J2000), with a depth of - 12mmag. 
The period of the signal was P = 4.1878 days, while 
the relative duration (first to last contact) of the tran- 
sit events was q « 0.027, which is equivalent to a total 
duration of Pq « 2.6 hours (see Fig. la). 

3. FOLLOW-UP OBSERVATIONS 

3.1. Reconnaissance Spectroscopy 

In order to exclude the possibility of a false planetary 
detection, due to the misinterpretation of a transit-like 
signal caused by another astrophysical scenario (such 
as an F + M dwarf system) , we observed the candi- 
date HTR176-002 with the CfA Digital Speedometer 
(Latham 1992) on the FLWO 1.5 m Tillinghast reflec- 
tor. We acquired four spectra between 2007 January and 
March, each with an individual precision of 0.5 km s^^. 
The observations showed a mean radial velocity of 7 = 
— 10.6kms~^ with an rms of 0.3kms~^, therefore rul- 
ing out a low-mass stellar companion (but not a triple 
system) , which would cause significantly higher RV vari- 
ations. The spectroscopy also yielded an estimate for the 
projected rotational velocity and surface gravity of the 
star. 

3.2. High S/N Spectroscopy and Subsequent Analysis 

We obtained high resolution and high signal-to-noise 
spectra with the Keck-I telescope and HIRES instrument 
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Fig. 1.— (a:) The light curve of XO-5 with all 14290 points 
taken in the I band, by the telescopes HAT-6, HAT-7, HAT-9 
and WHAT. The light curve is folded with the period of P = 
4.187757 ±0.000011 days (which is the result of the fit described in 
§ 5). The superimposed curve shows the best fit model, neglect- 
ing limb darkening, (b:) Unbinned instrumental Sloan z band 
and i band follow-up transit photometry light curves acquired 
with KeplerCam on the FLWO 1.2 m telescope on 2008 January 
2 (Mr = 0, z band), January 19 (A^tr = 4, 2 band), February 9 
(ATtr = 9, i band) and March 26 (JVtr = 20, i band). Superimposed 
are our best-fit transit models (§ 5). 

(Vogt et al. 1994). We acquired 17 exposures with the io- 
dine cell, and an additional iodine-free "template" . The 
measurements were made between 2007 March 27 and 
2008 May 17. The purpose of these observations was 
threefold: i) to obtain high precision radial velocity (RV) 
measurements, ii) to characterize the stellar properties, 
and iii) to check for spectral line bisector variations as 
an indication of blends. These steps are discussed in the 
following paragraphs. 

As regards measuring the RV variations, the super- 
imposed dense forest of I2 absorption lines enables us 
to obtain an accurate wavelength shift compared to the 
template observation (Marcy & Butler 1992; Butler et al. 
1996). The final RV measurements and their errors are 
listed in Table 1. The folded data, with our best fit (see 
§ 5) superimposed, arc plotted in Fig. 3, upper panel. 

The stellar atmosphere parameters were dcitermincd 
using the iodine-free template spectrum. The spectral 
modeling was performed using the SME software (Valenti 
& Piskunov 1996), with wavelength ranges and atomic 
hue data as described by Valenti & Fischer (2005). We 
obtained the following initial values: effective tempera- 
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ture Toff = 5505 ± 70 K, surface gravity logg^ — 4.61 ± 
0.10 (cgs), iron abundance [Fe/H] = +0.16 ±0.06, and 
projected rotational velocity z;sini = 0.7± 0.5kms~^. 

3.3. Photometric follow-up observations 

We obtained follow-up photometric observations 
on four nights using the KeplerCam CCD on the 
FLWO 1.2 m telescope through Sloan z and i bands. The 
observations were performed on 2008 January 2 (partial 
transit), January 19 (full transit), February 9 (partial) 
and March 26 (full), with the total number of object 
frames being 114, 428, 268 and 521, respectively. The 
integration times used at these nights were 45, 30, 30 
and 15 seconds, respectively while the readout and stor- 
age required an additional ^12 seconds per frame. The 
typical rms of the follow-up light curves was 2 mmag at 
the above cadence. 

We performed aperture photometry on the calibrated 
frames, using an aperture series that ensures optimal 
flux extraction. Details on the astrometry, photome- 
try, decorrelation for trends, etc., have been discussed 
in, e.g., Bakos et al. (2007). The light curves are plot- 
ted in the lower panel of Fig. 1, superimposed with the 
best-fit transit light curve model (see § 5). 

4. BLEND ANALYSIS 

A stellar eclipsing binary that is unresolved from a 
bright source would manifest itself as a blended system 
with shallow photometric transits, and with RV varia- 
tions that are of the same order of magnitude as one can 
expect from a planetary system (e.g. Queloz et al. 2001). 
We investigated whether such a blend is a feasible phys- 
ical model for HTR176-002 in two ways: by examining 
the spectral line bisectors, and with a detailed modeling 
of the light curve under various possible blend scenarios. 

For a blended eclipsing binary, in addition to the de- 
crease in the observed RV amplitude, the spectral lines 
would be distorted, as quantified by the "bisector spans" 
(see Torres et al. 2005, 2007). If the bisector span varia- 
tions correlate with the orbital phase, or the magnitude 
of these variations is comparable with the RV amplitude, 
then the system is likely to be a false positive (hierarchi- 
cal triple or chance alignment with a background binary) 
rather than a single star with a planetary companion. In 
order to rule this out we derived the bisector spans by 
cross-correlating the iodine-free ranges of the obtained 
spectra against a synthetic template spectrum. We found 
that the standard deviation of the bisector spans is ap- 
proximately ^ 60 m s^^ , which is comparable to the mag- 
nitude of the RV variation itself {K = 144.9 ± 2.0 ms~^; 
see Table 3). The large bisector variations discouraged us 
from publication even after the first full transit follow-up 
light curve was obtained in January 2008, and we con- 
tinued acquiring high resolution spectroscopy to establish 
whether there is any significant correlation between the 
bisectors and the orbital phase (or equivalently, with the 
actual RV values). In Fig. 4 we display our measure- 
ments of the bisector spans as the function of both the 
RV and the RV residuals from the best fit"'^^. There is no 

As wc will discuss later, our finally accepted best fit values 
were derived by including a decorrelation factor against this bisec- 
tor span correlation. In the plot the RV residuals are shown before 
subtracting this correlation term. 



statistically significant correlation between the velocities 
and bisector variations, as would be expected for a blend. 
However, there is apparently a correlation between the 
RV residuals and the bisector spans. This could be due to 
activity on the star (e.g., spottedness), where the activ- 
ity (if periodic) causes both RV and bisector variations, 
but in a way that is not commensurate with the orbital 
period of the companion. We exploit this correlation in 
the joint analysis of the RV and photometric data (see 
§ 5) where we show that the unbiased residual of the RV 
signal can be significantly decreased with the inclusion 
of an additional term proportional to the bisector spans. 

In order to rule out or confirm the importance of the 
stellar activity, we computed the Ca II emission index S 
(Noyes et al. 1984). The derived indices are also shown 
in Table 1. We found that the mean value of 5 = 0.16 ± 
0.02 is moderately low, and the correlations between the 
values of S and the radial velocity data or RV fit residuals 
are negligible (see also § 5.1). 

As a further way of assessing the true nature of the can- 
didate, we investigated possible blend configurations by 
performing light curve fits of our highest-quality follow- 
up photometry (data in the Sloan i band) following the 
procedures described by Torres et al. (2004). Briefly, we 
attempted to reproduce the observed photometric varia- 
tions with a model based on the EBOP binary-fltting 
program (Etzel 1981; Popper & Etzel 1981) in which 
three stars contribute light, two of which form an eclips- 
ing binary with the orbital period found for XO-5. The 
light from the third star (the candidate) then dilutes the 
otherwise deep eclipses of the binary, reducing them to 
the level observed for HTRl 76-002 (-1.2% depth). The 
properties of the main star were adopted from the results 
of our analysis below, and those of the binary compo- 
nents (mass, size, brightness) were constrained to satisfy 
representative model isochrones. We explored all pos- 
sible combinations for the binary components, and de- 
termined the best fits to the light curve in a chi-square 
sense. 

The case of a hierarchical triple (all stars at the same 
distance) yielded an excellent fit to the photometry (see 
top curve in Figure 2), but implies an eclipsing binary 
with a primary that is half as bright as IITR176-002 
itself. This is clearly ruled out by our Keck spectra 
and even our Digital Speedometer spectra, both of which 
would show obvious double lines. 

We then considered scenarios in which the eclipsing bi- 
nary is in the background (which would make it fainter) , 
and is spatially unresolved. Because the proper mo- 
tion of the candidate is relatively small (~30 mas yr~^; 
Monet et al. 2003), the chance alignment would remain 
very close for decades, precluding the direct detection of 
the binary in archival photographic images such as those 
available from the Digital Sky Survey. For convenience 
we parametrized how far behind the eclipsing binary is 
placed relative to the candidate in terms of the difference 
in distance modulus. Am, and we explored a wide range 
of values. As an example, we find that for Am = 4 (bi- 
nary about 1.7 kpc behind) the best fit yields a relative 
brightness for the binary of only 5%, which is at or below 
our detection threshold of 5-10% from the Keck spectra. 
However, the ingress and egress are clearly too long given 
the quality of our photometry (Figure 2, bottom curve). 
For a smaller separation of Am = 2 (binary some 500 pc 
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Fig. 2. — Blend modeling for XO-5, based on our Sloan i-band 
photometry. As examples we show the best fits corresponding to 
three different blend scenarios, with the bottom two displax;ed ver- 
tically for clarity. Top: Model corresponding to a hierarchical 
triple (see text), which is ruled out because the implied brightness 
of the eclipsing binary is so large (~50%) that our spectra would be 
double-lined. Middle: Model corresponding to a chance alignment 
with a baekgrouud eclipsing binary, in which the distance modulus 
difference between the binary and the candidate is Am = 2. The 
ingress and egress are already seen to be too long, and the fit im- 
plies a relative brightness of --^20% that would be easily detectable 
spectroscopically. Bottom: Chance alignment model with Am = 4 
in which the binary is much fainter (~5%), but the best-fit model 
does not match the observations well. These simulations rule out 
background blend scenarios. 

behind) the fit is somewhat better, though still visibly in 
disagreement with the observations (Figure 2, middle), 
and the relative brightness increases to 20%, which we 
would have noticed. Additional tests changing the incli- 
nation angle from the edge-on configurations considered 
above to lower angles did not alleviate the discrepancies. 

The above modeling rules out both a hierarchical triple 
and a background eclipsing binary as possible alternate 
explanations for the photometric signals we detect. This, 
combined with the lack of any clear correlation between 
the bisector spans and the radial velocities, constitutes 
compelling evidence of the planetary nature of HTR176- 
002 = XO-5, and convinces us that the scatter in the 
bisector spans described above is intrinsic to the star. 

5. ANALYSIS 

In this section we describe briefly our analysis yielding 
the orbital, planetary and stellar parameters for the XO- 
5 system. 

5.1. Light curve and radial velocity analysis 

For the initial characterization of the spectroscopic or- 
bit, we fitted a Keplerian model to the Keck RV data, 
allowing for eccentricity by including as adjustable pa- 
rameters the Lagrangian orbital elements k = ecosw 



and h = esintn, in addition to a velocity offset 7, the 
semi-amplitude K and the epoch E. The period P was 
held fixed at the value found from the HATNet light 
curve analysis (from BLS, see above). We found that k 
and h are insignificant compared to their uncertainties 
(k = -0.003 ± 0.029, h = -0.009 ± 0.023), suggesting 
that the orbit is circular However, in the determination 
of the orbital and stellar parameters, we incorporated the 
uncertainties yielded by the k and h orbital elements. 

We proceeded next with a joint fit using all data sets, 
namely, the HATNet discovery light curve, the FLWO 
1.2 m follow-up light curves, and the Keck radial ve- 
locities along with the initial estimates of the spectro- 
scopic properties derived through the SME analysis. The 
follow-up light curves were modeled using the analytic 
formalism of Mandel & Agol (2002), assuming quadratic 
limb darkening. The limb darkening coefficients 71.2, 
7i,i) 72, z and 72. i were taken from Claret (2004), in- 
terpolating to the values provided by the initial stellar 
atmospheric analysis § 3.2. The adjusted parameters for 
the joint fit were Tc^-27a, the time of first transit center 
in the HATNet campaign, Tc^20^ the time of the tran- 
sit center at the last follow-up (on 2008 March 26), m, 
the out-of-transit magnitude of the HATNet light curve 
in the I band, the semi-amplitude of the radial velocity 
K, the velocity offset 7, the Lagrangian orbital elements 
k and h, the fractional planetary radius p = i?p/i?^, 
the square of the impact parameter 6^, the quantity 
C/R, = (27r/P)(a/ii.)(l - 62)-i/2,/r3^(i + _ 
which is related to the duration of the transit-'^^ as 
2dur = 2(C/i?j,)~^, and the out-of-transit magnitudes 
Wc.o, 'Tic. 4, "T-c.g and mc,2o for the four follow-up light 
curves. See Pal et al. (2008) for a detailed discussion 
about the advantages of this set of parameters. The ini- 
tial values were based on the BLS analysis, and our ini- 
tial characterization of the orbit. To obtain the best-fit 
values, we utilized the downhill simplex algorithm (se(^ 
Press et al. 1992). The uncertainties and the correla- 
tions were determined using the Markov Chain Monte- 
Carlo method (Ford 2006) which yields the a posteriori 
distribution of the adjusted values. 

As mentioned in § 4, we found that there is a significant 
correlation between the RV residuals and the bisector 
spans. This suggests it might be possible to improve the 
RV fit by including an additional term to account for 
this correlation. We therefore expanded the model for 
the velocity variation to 

Vi=j + K.R\o(^ ^""^^'p ,k,h^+ CBsbi (1) 

where RVo(-, •, •) represents the base function for the ra- 
dial velocity variations'^ and 6, is the actual bisector 
span variation for the i-th measurement. We found that 
when omitting the last term the unbiased residual is 

8.8 ms^', whereas its inclusion leads to decreased residu- 
als of 4.6 ms^', nearly a factor of two better. We tested 

Here duration is not the total duration between the first and 
last contact but defined as the interval between the instances when 
the center of the planet crosses the the limb of the stars inward and 
outward. 

This function has three arguments: the mean longitude 
measured from the transit center and the two Lagrangian or- 
bital elements k and h. It is easy to show that if fc = /i = 0, 
RVo(A,0,0) = -sin(A). 
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TABLE 1 

Relative radial velocity, bisector span and stellar 

ACTIMT^' f.SblEASrREirENT^ OF XO-o 
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0.1530 
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226.59 
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2.79 


20.77 


42.04 


0.1589 


216.76639 


294.36 


2.70 


81.43 


20.90 


0.1549 




0.00 


3.37 


—30.10 


45.47 


0.1876 


248.77938 


76.84 


3.83 


-34.80 


48.41 


0.1855 


249.78531 


268.26 


3.36 


12.13 


36.59 


0.1548 


251.78153 


5.41 
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0.1459 


428.02826 


8.16 


3.02 
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94.90 
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3.41 
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224.82 


7.25 
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0.1577 


548.81658 


79.96 


3.08 
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0.1595 


548.89652 


65.98 
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24.33 


0.1576 


602.74168 


193.74 


2.61 


59.89 


25.39 


0.1578 
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2.83 


60.49 


25.44 
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''This spectrum turned out to be severely contaminated by moon- 
light; however, the corresponding RV is unaffected. 
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Fig. 3. — (a) Radial-velocity measurements from Keck for XO-5, 
along with our orbital fit (sec § 5). The centcr-of-mass velocity 
7 and the correlation correction for the bisector span variations 
has been subtracted, (b) Phased residuals after subtracting the 
orbital fit (also sec § 5). The rms variation of the residuals is 
about 4.6 ms^-'-. (c) Bisector spans (BS) for 16 of the Keck spectra 
(including the iodine- free template). Note that the scales of the 
panels are the same. 

also whether the inclusion of a similar term in Eq. (1) 
proportional to the stellar activity index (with a coef- 
ficient Cs- index) provides any further improvement in 
the fit, but found that it actually degrades the residu- 
als slightly. The final orbital and planetary parameters 
(and their uncertainties) derived in this paper are based 
on the above discussed radial velocity model function 
decorrelated against the bisector variations. 

5.2. Stellar and planetary parameters 

The stellar parameters were determined in an iterative 
way as follows. As pointed out by Sozzetti et al. (2007), 
the stellar density is a better luminosity indicator than 



TABLE 2 
Stellar parameters for XO-5 

Parameter Value Source 

Tgff (K) 5370 ± 70 SME* 

[Fe/H] +0.05 ± 0.06 SME 

i;sini (kms-i) 0.7 + 0.5 SME 

M* (Mq) 0.88 + 0.03 Y^+LC+SME^ 

Ri, {Rq) 1.08 + 0.04 Y^+LC+SME 

log (cgs) . . . 4.31 + 0.03 Y^+LC+SME 

(Lq) 0.88 + 0.09 Y2+LC+SME 

Mv (mag) .... 5.06 + 0.12 Y^+LC+SME 

Age (Gyr) .... 14.8 + 2.0 Y^+LC+SME 

Distance (pc) . 260 + 12 Y^+LC+SME 

^SME = 'Spectroscopy Made Easy' package 
for analysis of high-resolution spectra Valenti & 
Piskunov (1996). Sec text. 

''Y^+LC+SME = Yalc-Yonsci isochroncs (Yi et al. 
2001), light curve parameters, and SME results. 
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Fig. 4. — Bisector span variations as a function of the RV (left 
panel) and RV fit residual (right panel). The rigfit panel sfiows the 
fit residuals when the correlation term was not included in the fit. 
Note that on the graphs the horizontal scales are not the same. 



the spectroscopic value of log g^.. In a first order approx- 
imation the density is related to the observable quanti- 
ties P and a/Ri, as = (37r)G~^P~^(a/i?^)^. We used 
the values of T^s and [Fe/H] from the SME analysis, to- 
gether with the distribution of (derived from a/Ri,) 
to estimate the stellar parameters from the Yonsei-Yale 
evolution models, as published by Yi ct al. (2001) and 
Demarque ct al. (2004). This resulted in a posteriori 
distributions of those stellar parameters, including the 
mass, radius, age, luminosity and colors. Prom the mass 
and radius distributions, we obtained a new value and 
uncertainty for the stellar surface gravity: 4.31 ± 0.03. 
Since this value is significantly smaller than the previ- 
ous value based on the SME analysis § 3.2, we repeated 
the atmospheric modeling by fixing the surface gravity 
to the new value (4.31 ± 0.03), and allowing only the 
metallicity and effective temperature to vary. This next 
iteration of the SME analysis yielded T^s = 5370 ± 70K 
and [Fe/H] = +0.05 ± 0.06. Based on these new at- 
mospheric parameters, the limb darkening coefficients 
were re-calculated and we repeated the joint fit for the 
light curve and RV parameters, followed by the stellar 
evolution modeling once; again, in the same way as dis- 
cussed earlier. In this iteration the surface gravity barely 
changed (log 5* = 4.33 ±0.04), so the stellar previous pa- 
rameters were accepted as final (Table 2). In Fig. 5, we 
plot the evolutionary isochrones as the function of the 
effective temperature and both the stellar surface grav- 
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Fig. 5. — Stellar evolution isocliroiics from the Yoiisei-Yalc mod- 
els, corresponding to ages between 2 and 14 Gyr (in steps of 1 Gyr), 
as a funetion of both surfaee gravity (top) and normalized semima- 
jor axis a/R* (bottom). In the top panel the isoehrone metallieity 
([Fe/H] = +0.16), spectroseopic surfaee gravity, and temperature 
are from our initial SME analysis, the latter two shown with l-cr 
and 2-iT confidence ellipsoids. In the lower panel the metallieity 
([Fe/H] = +0.05), temperature, and a/R* are from the iterative 
analysis described in the text. Note that the latter quantities result 
in a significantly different evolutionary state for the star. 

ity and a/il* (these are used as luminosity indicators). 
The temperature, surface gravity and relative semimajor 
axis values discussed here are also superimposed on these 
isoehrone plots. 

The results from this second global fit to all the avail- 
able data (photometry, radial velocities) are listed in 
Table 3. In addition, values for some auxiliary pa- 
rameters in this fit are: -270 = 2453338.22311 ± 
0.00236 (BJD), Tc.ao = 2454552.67174 ± 0.00029 (B.JD), 
m = 11.33042 zbO.OOOlOmag and the Keck velocity offset 
is 7 = 0.8±0.1ms~^. The best-fit values and uncer- 
tainties for the fitted parameters are straightforward to 
obtain from the MC distributions. These, in turn, lead 
to the planetary parameters and their uncertainties by 
using a direct combination of the a posteriori parameter 
distributions of the light curve, radial velocity and stellar 
parameters. We find that the mass of the planet is Mp = 
1.059 ± 0.028 Mj, the radius is Rp = 1.109 ± 0.050 i?j 
and its density is Pp = 0.96toliigcm~^. These quanti- 



ties are also collected in Table 3. The correlation coeffi- 
cient C'{Mp, Rp) between the planetary mass and radius 
is listed as well. We also estimated the individual transit 
centers of the four follow-up light curves, by adjusting 
only the the light curve parameters (iip/i?*, b"^, C/-R*) 
out-of-transit magnitudes) while the transit centers were 
not constrained by a given epoch and period. We ob- 
tained that the individual transit centers do not differ 
significantly from the interpolated transit centers (de- 
rived from the results of the joint fit), i.e. the available 
data do not show any signs for transit timing variations. 
The independently fitted transit centers for the events 
= 0, Ntr = 4 and Ntr = 20 differ from the hnearly 
interpolated values by less than 1.5-0-, and the difference 
at the event A'tr = 9 is nearly 2.3-a. The independently 
fitted and the interpolated transit instants are shown in 
Table 4. 

Using our best fit model, we also checked the ampli- 
tude of the out-of-transit variations of the HATNet light 
curve, by performing a Fourier analysis on the fit resid- 
uals. We found no significant variation in the stellar 
flux, and all Fourier amplitudes were less than O.Tmmag. 
This estimation gives an upper limit for the stellar activ- 
ity, and is in line with the small S values derived from 
spectroscopy (S < 0.186, see Table 1). It is somewhat 
surprising that in spite of the small activity based on the 
spectroscopic S index, the light curve out-of-transit vari- 
ation, and the low vsini rotational velocity of the star, 
the bisector spans exhibit such a large scatter. 

The Yonsei-Yale evolutionary models also provide the 
absolute magnitudes and colors for different photomet- 
ric bands. We compared the V — I model color with 
the observed TASS color (see Droege et al. 2006). Since 
{V-I)yy = 0.815±0.016and (V-/)tass = 0.82±0.09, 
we conclude that the star is not significantly affected 
by interstellar reddening (also note the Galactic lati- 
tude of XO-5, which is b = 26°. 9). Therefore, for 
the distance determination we use the distance modu- 
lus Vtass — My = 7.18 ± 0.13, which corresponds to 
d = 260±12pc. 

6. DISCUSSION 

In this paper we have described our independent de- 
tection of the transiting planet X0-5b using the HAT- 
Net observations. A significant component of our effort 
has been to examine possible astrophysical false posi- 
tives and to model the data in detail in order to rule 
them out. In this way we have provided new and cru- 
cial support for the planetary nature of the object. We 
also present refined values for the system parameters. It 
is reassuring that the planetary parameters in B08 and 
this work are consistent within l-cr. This, however, is 
somewhat coincidental, since the stellar parameters are 
quite different. Based on our SME analysis, we derive a 
lower effective temperature {Tes= 5370 ± 70K as com- 
pared to 5510 ± 44K in B08), and a lower metallieity 
([Fe/H] = +0.05 ± 0.06 vs. 0.25 ± 0.03). The difference 
is attributed to our iterations on the SME analysis and 
the transit-fitting; using the n/i?* based mean stellar 
density as a luminosity indicator, and fixing the corre- 
sponding logs'* ill the SME analysis (i.e. solving only 
for [Fe/H] and Teg). We derive a smaller stellar mass: 
0.88 ± 0.03 Mq vs. 1.0 ± 0.03, M©, based on the same Yi 
et al. (2001) isochrones. Due to the high precision photo- 
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TABLE 3 
Orbital and planetary parameters 



Parameter 



Value 



Fig. 6. — (a): Safronov number vs. equilibrium temperature for the known transiting extrasolar planets. XO-5b is marked with a larger 
dot and it is located at the upper envelope of the Class I distribution of planets (XO-5b > 0.05). (b): Surface gravities as the function of 
the orbital period for the known extrasolaj: planets. With its relatively high surface gravity and orbital period, XO-5b falls slightly off the 
main distribution. 

metric and RV data, we are able to refine the planetary 
and orbital parameters of the system, and decrease the 
uncertainties typically by a factor of ~ 2 — 3. 

Based on the models of Liu ct al. (2008), after rc- 
scaling the semi-major axis to match the insolation flux 
X0-5b would have if it orbited a G2V dwarf {uequiv = 
0.05313 AU), the measured mass and radius of X0-5b re- 
quire a small core to be consistent with theory even if no 
internal heating is assumed. Using the work of Fortney et 
al. (2007), X0-5b is consistent with a 300 Myr old planet 
with a 50 core, a 1 Gyr old planet with a 25 core, 
or a 4.5 Gyr planet with a core smaller than 10 mass. 
The incident flux on X0-5b is ~ 4.83 • 10^ergs~^cm~^. 
This corresponds to a pL class planet, based on the defl- 
nitions of Fortney et al. (2008), although it falls fairly 
close to the transition area between the pL and pM 
classes. 

We confirm that the planet has a remarkably high 



Light curve parameters 

P (days) 4.187757 ±0.000011 

E (BJD) 2454552.67168 ± 0.00029 

Ti4 (days)-'' 0.1307 ± 0.0013 

Ti2 = T:i4 (days)'' 0.0175 ± 0.0013 

C/R* (day-i) 17.779 ± 0.091 

a/Ri, 9.67 ±0.35 

Rp/Ri, 0.1050 ±0.0009 

b = acosi/Ri, 0.562too5i 

i (deg) 86.7 ± 0.4 

Spectroscopic parameters 

K (ms-i) 144.9 ±2.0 

Cbs 0.125 ±0.025 

Cs-index (adopted) 

fc s e cos a; +0.008 ± 0.010 

/isesino; +0.010 ± 0.013 

Planetary parameters 

Mp (Mj) 1.059 + 0.028 

Rp (Rj) 1.109 ±0.050 

C{Mp,Rp) 0.23 

Pp (gcm-3) 0.96l° it 

a (AU) 0.0488 ± 0.0006 

log gp (cgs) 3.33 ± 0.04 

Teq (K) 1221 ± 27 

e 0.105 ±0.005 

T14: total transit duration, time between first and last 
contact; T12 = T34: ingress/egress time, time between first 
and second, or third and fourth contact. 



TABLE 4 

Individual transit center measurements 



Event 


TciBJDY 


Tc{BJDf 


# 


(2,454,000+) 


(2,454,000+) 





468.91868 ± 0.00181 


468.91666 ± 0.00028 


4 


485.66932 ± 0.00058 


485.66768 ± 0.00028 


9 


506.60475 ± 0.00057 


506.60645 ± 0.00027 


20 


552.67152 ± 0.00041 


552.67174 ± 0.00029 



"Derived frp, the individually fitted the transit cen- 
ters while the other light curve parameters were con- 
strained to be equal. 

''Derived by interpolation from the joint fit results, 
assuming a constant period. 



Safronov number, 6 = l/2{VesclVorbf = 0.105 ± 0.005, 
placing it at the high end of the Class I planets as defined 
by Hansen & Barman (2007). The plot of the Safronov 
numbers for the known TEPs as a function of equilibrium 
temperature is displayed on Fig. 6a. We also confirm that 
X0-5b has an anomalously high surface gravity, as com- 
pared to other TEPs with similar period (Southworth, 
Wheatley & Sams 2007). 

Altogether, XO-5 appears to be an interesting system 
exhibiting a number of anomalies including non-trivial 
bisector span variations, and anomalously high Safronov 
number and surface gravity. Future observations and 
theoretical work are required to understand these prop- 
erties. 
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